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Abstract 

This paper deals with Poisson processes on an arbitrary measurable space. Using a 
direct approach, we derive formulae for moments and cumulants of a vector of mul- 
tiple Wiener-Ito integrals with respect to the compensated Poisson process. Second, 
a multivariate central limit theorem is shown for a vector whose components ad- 
mit a finite chaos expansion of the type of a Poisson U-statistic. The approach is 
based on recent results of Peccati et al. combining Malliavin calculus and Stein's 
method, and also yields Berry-Esseen type bounds. As applications, moment for- 
mulae and central limit theorems for general geometric functionals of intersection 
processes associated with a stationary Poisson process of £;-dimensional flats in K 
are discussed. 

Key words: Berry-Esseen type bounds; central limit theorem; intersection process; multi- 
ple Wiener-Ito integral; Poisson process; Poisson flat process; product formula; stochastic 
geometry; Wiener-Ito chaos expansion 

MSC (2010) subject classifications: Primary: 60D05; 60H07 Secondary: 60F05; 60G55 

1 Introduction 

Throughout this paper we denote by (X, X) a measurable space, equipped with a er-finite 
measure A ^ 0. A classical result by Wiener and Ito says that if F = F(rj) is a square 
integrable function of a Poisson process rj on (X, X) with intensity measure A (cf. [TJ 
Chapter 12]), then F can be represented as an orthogonal L 2 -series 

oo 

F = EF + Y,In(fn), (1-1) 

n=l 



*Institut fiir Stochastik, Karlsruhe Institute of Technology, 76128 Karlsruhe, Germany. E-mail: 
gucnter.last@kit.edu 

^Department of Mathematical Sciences, University of Bath, Bath BA2 7AY, United Kingdom. E-mail: 
m.d.penrose@bath.ac.uk 

■flnstitut fiir Stochastik, Karlsruhe Institute of Technology, 76128 Karlsruhe, Germany. E-mail: 
matt hias . schulte@kit.edu 

^Department of Mathematics, University of Osnabriick, 49076 Osnabruck, Germany Current 

address: Faculty of Mathematics, Ruhr University Bochum, 44801 Bochum, Germany E-mail: 
christoph.thaclc@rub.de. 



1 



where I n (f n ) is the multiple (n-fold) Wiener-Ito integral of a certain symmetric function 
f n : X n — > K that is square integrable with respect to A n . For a proof of this result without 
any further assumptions on the measure space (X, A", A) we refer the reader to [TUJ. It 
turned out that the chaos expansion (11. ip is useful for many purposes. For instance, it 
serves as mathematical basis for Malliavin calculus of variations on the Poisson space and 
can be used to formulate and to prove central limit theorems, see [T2l [T4] . 

The present paper deals with multivariate Poisson functionals (ij , • • • , Ft), ^ G N, 
where each component is of the form F^ = F^(r] t ) and 7] t is a Poisson point process 
with an intensity measure of the form X t = tX. Here, we are interested in the asymptotic 
regime that arises when the intensity parameter t tends to infinity. Under the additional 
assumption that each f} is a U-statistic of the Poisson process rjt, we prove formulae for 
the joint moments and cumulants and a multivariate central limit theorem. 

The assumption that the functionals F^ are Poisson U-statistics implies that their 
chaos expansions ( II. ip terminate after finite numbers of terms, which is convenient for the 
application of Malliavin calculus. Univariate central limit theorems with bounds on the 
Wasserstein distance for Poisson functionals with finite Wiener-Ito chaos expansions and, 
in particular, Poisson U-statistics are derived in [H|9l[T6] using a general result from [12]. 
Our multivariate counterpart rests on a multivariate analogue in [TJ]. By using a new 
truncation argument and the special form of the Poisson functionals, our approach avoids 
technical computations involving the product formula for multiple Wiener-Ito integrals 
that is used in the works mentioned previously. 

As an application we study in detail geometric functionals of the intersection process of 
order m 6 {1, . . . , d} of a stationary Poisson fc-flat process in M. d . We thereby considerably 
extend the results available in the literature (H |5] for the number of intersections and the 
intersection volume. In our theory we can allow for very general geometric functionals, 
for example, we do not require them to be additive, translation-invariant or homogeneous. 
Furthermore, our central limit theorems are quantitative in the sense that they provide 
rates of convergence (with respect to a suitable distance). 

As well as in the context described above, functionals of the type we consider show 
up in a number of other problems arising in stochastic geometry. Examples to which 
our theory may directly be applied include subgraph counting functionals for random 
geometric graphs [HJ [9j [15], the Euler characteristic of random simplicial complexes [3], 
the number of fc-clusters [TJ, random graph statistics of high-density Boolean models [9], 
and proximity functionals for non-intersecting /c-flat processes |18j . 

Another direction this paper deals with is formulae for mixed moments and cumu- 
lants, which in turn are based on identities for mixed moments and cumulants of multiple 
Wiener-Ito integrals. We develop a quick approach to prove these formulae that have pre- 
viously appeared in different generality in [T3| [19] . The novelty of our proof is that it only 
makes use of elementary properties of the Poisson process (mainly the multivariate Mecke 
formula) and some combinatorial arguments, and deals directly with the expectation. In 
this way it avoids requiring the involved chaos expansion of such products. 

The text is structured as follows. In Section [2] we collect some basic definitions and 
background material. The moment formulae are presented in Section [3] while Section 
H] deals with their asymptotic behaviour. Our multivariate central limit theorems for U- 
statistics of the underlying Poisson process are the content of Section [5] whereas in Section 
E]our results are applied to Poisson fc-flat processes in M d . 
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2 Preliminaries 



In this paper all random objects are defined on a probability space (Q, J 7 , P). We interpret 
the Poisson process rj as a random element in the space N := N(X) of integer- valued 
(including +00) a-fmite measures /1011X equipped with the smallest u-field Af making 
the mappings fi t— > [i(B) measurable for all B G X. For m G N and \i G N we define a 
measure on X m by 



m— 1 m— 2 



• • • / i B (gi, . . . ,x m )(ji - 5 Xj j (dx m ) f/x - <j Xi j (dx m _i) 
J 3=1 3=1 

. . . (fi - S xi )(dx 2 )fi(dx 1 ), 

where S x is the Dirac measure located at a point x G X. If \x is given as /x = ^ 7 for 
some countable index set / and Xj G X, % G /, then 



/ /d/i( m )= ^ /(.r, : r„J. 



h,...,i m €l 



where / is any non-negative measurable function on X' m and where the superscript 7^ 
indicates that we sum over m-tuples of disjoint indices. 

We will use the multivariate Mecke-formula (see e.g. [TO] ) 



E / h(r],xi, . . . ,x m )r/ (m) (d(xi, . . . ,x m )) 

J f (2.1) 

= E / h(r) + 5 xi + ... + 5 Xm ,x u ...,x m ) A m (d(xi, . . .,x m )) t 

which holds for all h : N x X m — > K for which one (and then also the other) side makes 
sense. 

For any integer n > 1 let Lj,(A n ) denote the set of all measurable and symmetric 
functions / : X n — > K that are integrable with respect to A n . For / G L l s (X n ) define the 
pathwise multiple Wiener-Ito integral by 

W) ■= E^ 1 )" 1 ' 1 [[ f(xu---,Xn)v m (dxj)\ n - lJ Kdxj«), (2.2) 

JC[n] JJ 

where [n] :—{!,..., n}, J c := [n] \J, xj := (xj)j e j and where \ J\ denotes cardinality of J 
(the inner integral is interpreted as f(xi, . . . , x n ) in the case where J = 0). By (12. ip . I n (f) 
is a well-defined integrable random variable with EJ n (/) = 0. If / G L](\ n ) n L 2 (A n ), 
this pathwise definition coincides with the (classical) definition of the multiple Wiener- 
Ito integral for square integrable functions. (This fact can be derived, for instance, from 
equation (3.1) in [ID].) The multiple Wiener-Ito integral I n {f) for (symmetric) / G L 2 (X n ) 
is defined by an extension of the pathwise definition for L 1 -functions to the space of all 
square integrable random variables. It has also mean zero and satisfies the orthogonality 
and isometry relations 

EI m (g)I n (h) = l{m = n} n\ (g, h) n , m,n>l, (2.3) 
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for all (symmetric) g £ L 2 (A m ) and h £ L 2 (X n ), where ( • , ■ ) n denotes the scalar product 
in L 2 (A n ). 

For x £ X the difference operator D x is given as follows. For any measurable F : N — y 
R the function D X F on N is defined by D x F(ji) := F(fi + 5 X ) - F(ji), /i £ N. Forn > 2 
and (xi, . . . , x n ) £ X n we define a function .D™ Xn F : N — > R by an iterated application 
of the difference operator D, that is, inductively by 

£> n F := D 1 D n ~ l F 

where .D* := -D x and Z?°.F := F. Under the assumption E,F(n) 2 < oo it was proved in 
[TU] that Z?£ x„F(v) is integrable for A n -a.e. (x l7 . . . , x n ) and that 

T n F(x 1 , . . . , £ n ) :=ED« ( ....^Ffo), . . . , x n ) £ X", 

defines a symmetric function in L 2 (A n ). Moreover, we have the Wiener- Ito chaos expansion 

F(r / ) = EF + ^-/ n (T n F), (2.4) 

n=l 

where the series converges in L 2 (P). Hence holds with / n = -^T n F. 

3 Moments and cumulants 

Let n £ N. A subpartition of [n] is a family of disjoint and non-empty subsets of [n], which 
we call blocks. A partition of [n] is a subpartition a of [n] such that Uj gcr J = [n]. We 
denote by Il n (respectively LT*) the system of all partitions (respectively subpartitions) 
of [n\. The cardinality of a £ II* (i.e. the number of blocks of a) is denoted by \o~\, while 
the cardinality of Uj ecr J is denoted by \\o~\\. For any function / : X n — y M and a £ II* 
we define f a : X' CT ' +n ~" CT " — y M. by identifying the arguments belonging to the same J £ a. 
(The arguments x±, . . . , x\ a \ +n -\\ a \\ have to be inserted in the order of occurrence.) In the 
case n = 4 and a = {{2, 3}, {4}}, for instance, we have f a (xi, x%, £3) = f(%i, X2, %2, ^3)- 
Consider ^, m, . . . , ni £ N. Define n := ri\ + • • • + and 

Jj := {j : ni H hnj_i<j<niH hn^}, i = l,...,£. (3.1) 

Let 7r := { Ji : 1 < i < £} and let Il(ni, . . . , ni) C Il n (respectively II*(ni, . . . , ni) C II*) 
denote the set of all a £ Il n (respectively er £ II*) with | J D J'\ < 1 for all J £ 7r 
and all J' £ er. Let n> 2 (ni, . . . ,ni) (respectively n =2 (ni, . . . ,ni)) denote the set of all 
er £ Tl(ni, . . . ,ni) with \ J\ > 2 (respectively |J| = 2) for all J £ er. It is instructive 
to visualize the pair (ir,(x) as a diagram with rows Ji, . . . , Jg, where the elements in 
each J £ er are encircled by a closed curve; see [131 Chapter 4] for more details on such 
diagrams. Since the blocks of a er £ II(ni, . . . , ni) are not allowed to contain more than one 
entry from each row, the diagram (tt, er) is called non-flat in [13]. Any a £ II>2(ni, . . . , ni) 
induces a partition er* £ IL: er* is the finest partition of [£} such that two numbers i,j £ [£] 
are in the same block of er* if Ji and Jj are both intersected by the same block of er. Let 
n>2(rii; . . . ,ni) be the set of all er £ II>2(ni, . . . ,ni) such that |er*| = 1. 
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The tensor product (8>f =1 /t of functions : X n * — > R, i G {!,...,£}, is the function 
from X n to R which maps each (x±, . . . , x n ) to Yli=i m case that f 1 = ... = f i = f 

we write /® instead of ®f =1 /i- 

The jomt cumulant 7(^1, . . . , Xij of £ > 1 random variables Xi,...,Xt is defined as 



j(X 1 ,...,X e ) := (-0^ ^logE[exp[i(z 1 X 1 + --- + ^)] 



Z! = --- = Z e =0 



where i is the imaginary unit. This cumulant is well-defined if rijg/^i ^ s integrable for 
all / C [£}. The £-th cumulant of a single random variable X is defined by ji(X) := 
"f(X, . . . , X), where X appears £ times. 

The following result generalizes [TBI Corollary 7.4.1] and a consequence of [19] to 
the case of more general Poisson processes and integrands. In contrast to p31 [19], we 
allow that the intensity measure has atoms. Moreover, we avoid the assumption in [TBI 
Corollary 7.4.1] that the integrands are simple functions. While the results in [TBI [19] 
are derived via formulae for the Wiener-Ito chaos expansion of a product of multiple 
Wiener-Ito integrals, here we take a direct approach, which relies only on ( 12.1 p and some 
combinatorial arguments. 

Theorem 3.1. Let fa G L](\ ni ), i — 1, . . .,£, where £, n±, . . . , ri£ G N. Assume that 

dl/,|) CT dA |CT| < 00, a G n(ni, . . .,n £ ). (3.2) 
Then 



i 

Ej[l nt (f i )= Yl yCd/OadAH, (3.3) 

i=l <rGn>2(ni,...,n^) 

7(/» 1 (/i),...,/ B i(/^))= E /(Ci/i)adAW. (3.4) 



ff£n> 2 (ni,...,n ? ) 

Proof. We abbreviate / := ® i=1 fi- The definition (12. 2p and Fubini's theorem imply that 

n / «.(/i)=E(- i ) w " i/i / /(^.■•■^) 

i=i /cm 7 7 ( 3 - 5 ) 

r/' /nJl ')(d* /nJl ) . . .^^(da^A^'^), 



where 7 C := [n] \ I, and where we use definition (13.1 p of Jj. For fixed I C [n] we may 
split the above integration according to a G U*(ni, . . . , n<), where U,/ Gcr J = /. For any 
such cr we integrate (i.e. sum) over those (xi, . . . , x n ) satisfying £j = xj whenever i and j 
belong to the same block of a but not otherwise. By (12. ip applied with h(rj, yx, . . . , y\ a \) = 
f(x±, . . . ,x m ) taking Xi = yj for i in the j-th block of a, the contribution of a to the 
expectation of the right-hand side of (13. 5 p equals (— l) n ~H CT H J /„. dA'°"' +n ~" CT ". Therefore, 

- r 
Ej[l ni {f i )= Yl (-l) ft "W / / ff dAW+»-IWI. (3.6) 

i=i (Ten*(ni,...,%) 
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By assumption (13.21) all of these integrals are finite. Take a a £ U*(ni, . . . , ng) with | J| > 2 
for all J £ a and consider the set n 1 (a) of all 0\ £ U*(ni, . . . ,ng) such that a C o\ and 
| J| < 1 for all J £ (7i \ (7. (Note that a £ Ilita).) Observe that / dAl* 71 ^"' 71 " = 
/ / a dAH+'HMI for all ai £ IIi(ct). Moreover, for 

n ~ II " II > 1 we have that 

E (-ir lki11 = o. 

<rierii(<T) 

Since every r £ U*(ni, . . . , ni) has a unique c £ II*(ni, . . . , n$) with |J| > 2 for all 
J £ cr such that r £ ni(cr), we can partition II*(ni, . . . , nt) into the sets Ui(a), a £ 
Il*(ni, . . . ,ng) with |J| > 2 for all J £ cr. As shown above the sums over all a± £ Hi(cr) 
with \\a\\ < n vanish and only the integrals related to partitions a £ II>2(ni, . . . , ni) 
remain. Therefore, (13.61) implies the asserted identity (13. 3p . 

We now prove (|H3D by induction over £. Since j(l m {fi), ^(h)) = E/ ni (/i)/ n2 (/ 2 ) 
and n> 2 (ni,?7, 2 ) = n> 2 (ni, n 2 ), the identity (13.41) is true for £ = 2, see (I3.3p . For £ > 3, 
we obtain by the inversion formula expressing the £-th moment in terms of lower order 
cumulants (see e.g. [T3l Proposition 3.2.1]), formula (13.31) and the assumption of the 
induction 

l(l ni {fi),...,I nt {fa)) 

i 

=ElI^(/i)- E U'^'-Mmr') 

j=l 7rGn £ ,|7r|>l JStt (3- 1) 

= E ju^- E II E Ji^jfjUdx^. 

aGJI> 2 (ni,...,n e ) 7ren f ,|7r|>l Je?r aj&fl> 2 {J) 

Here I1> 2 (J) is defined in a similar way to n> 2 (rii, . . . Now we use the fact that 

every partition a £ Il> 2 (ni, . . . ,rii) determines (in the obvious way) uniquely a partition 
7r = cr* £ 11^ and a collection of partitions oj £ IT> 2 ( J), J £ ir, and vice versa. Combining 
this with Fubini's theorem, we have 

E /(n®^^ dAH 

c r £n>2(ni,...,rif),|(j* |>1 JGcr* 

E /^ dAH 

a"€n>2(ni,...,?i f ),|(T* |>1 

Hence in (13.71) only the partitions cr £ n> 2 (ni, . . . , ng) with \a*\ = 1 remain. In our 
notation, these are exactly the partitions in n> 2 (ni, . . . , ng). □ 

Remark 3.2. The assumption ( 13 .2p is obviously satisfied if is bounded and \ ni ({fi ^ 
0}) < oo for % — 1, . . . , £, which is the case for our examples in Section [6j But the assump- 
tion also holds under the following conditions. Consider, as in Theorem I3.1[ measurable 
functions /j : X n ' — > M, for i = 1, . . . ,£. Assume for any i that {/j ^ 0} C B n \ where 
B £ X satisfies X(B) < oo. For any i, if j i £ L e (\ n >) then f t £ L l (\ ni ) and flM} holds. 
The second assertion follows from the multivariate version of Holder's inequality. In fact, 



E II E j&XjMajdX™ = 
7ren £ ,|7r|>l JSvr o-j6fl>2(J) 
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if a G n(rii, . . . , rii), then 

y^LilfSadX^ <A(B)W~ ni y \f 1 \ e d\^-...-\(B)^ n ' J Iftl'dX*. 

Another sufficient condition for the assumptions of Theorem 13.11 is the existence of a 
function g G -^(A) PI L (X) such that < g® Ui for any i. In this case we have for 
cr G n(ni, . . . , Hi) that 

(CiI/.D^AH < / dA ..... / ^vi dA, 



where z'i, . . . , < £ are the cardinalities of the blocks of a. 

Example 3.3. Let / G L*(A 2 ) and consider Theorem 13.11 in the case I = 2, n\ = n 2 = 2 
and fi — f 2 — f. Then it is easy to see that assumption (13.21) requires / G L 2 (A 2 ) and 



|/(xi,x 2 )| A(dx 2 ) 



A(da;i) < oo. 



Formula (13.31) boils down to the isometry relation E/ 2 (/) 2 = 2(f,f) 2 - This shows that 
assumption (I3.2p is not necessary for (13.31) . 

Example 3.4. Let / G L>l(X 2 ), g G -^ 1 (A) and consider Theorem 13.11 in the case I = 3, 
n i = n 2 = 2, fi — J 2 — f, n 3 — 1 and f 3 = g. Assumption (13. 2p then requires / to satisfy 
the same integrability conditions as in Example 13.31 as well as 

(\f(x 1 ,x 2 )\ + /(xi,x 2 ) 2 )|5'(a:i)| A 2 (d(xi,x 2 )) < oo, 

\f(x 1 ,x 2 )f(x 2 ,x 3 )\(\g(x 1 )\ + \g(x 2 )\) X 3 (d(x 1 ,x 2 ,x 3 )) < oo. 

Formula (13.31) means that 

E[I 2 (f) 2 I 1 (g)]=4 / f(x 1 ,x 2 ) 2 g(x 1 )X 3 (d(x 1 ,x 2 )). 



Note that we do not need to assume g to be square-integrable with respect to A. 
Corollary 3.5. Let f n G Ll(X n ), n G N, and let i G N and assume that 

y"(®Lil/n,|) CT dA H < oo, cr g n(m,...,rii), nx,...,n t G N. 



(3.« 



Assume further that K(Y2^ =1 \I n {fn)\Y < oo. Then the l-th moment and the l-th cumulant 
of F := Er=i J «(/n) are 9 iven b V 



Ep£ = E E n®iiU,^\ (3.9) 

ni,...,n^gN crgn>2(ni,...,n f ) 

7^)= E E /(d/nJ.dAH. ( 3.io) 
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Proof. Let m G N and F m := h(fi) + ■ ■ ■ + I m {fm)- Expanding (F m ) 1 and using (13. 3p 
gives 

m 

EFi= Yl E(J ni (/ ni )....-J ni (/nJ) 

ni,...,n£=l 



E E /(Cl/nJ.dAH. 



rii,...,n f =l crGn> 2 (ni,...,n f ) 

Assumption E(J^L 1 |^n(/n)|) £ < oo and dominated convergence imply (13. 9p for the infi- 
nite case. By the multilinearity of joint cumulants and (13. 4p . 

m „ 

l£ (F m )= E E ]H=lfnX^- 

n 1 ,...,n l =l ff€ n> 2 (ni r .,Ti() 

Since 7^(-F m ) is a polynomial in the moments Ei 7 ^, j G {1, . . . ,£}, whose coefficients are 
independent of m (or F m ), we can again use dominated convergence to conclude the result 
(I3.f0p for the infinite case. □ 

Remark 3.6. If in Corollary 13.51 the number of non- vanishing functions /„ ^ is finite, 
then the assumption E(^^° =1 \I n (fn)\Y < oo is implied by ( 13. 8p . 

Let / G L 2 (X 2 ) be given by f(xi,x 2 ) := l{xi G B,x 2 G £?}, where A(£?) < oo. By 

J 2 (/) = n(5)(n(5) - 1) - 2n(B)A(B) + A(5) 2 . 

A straightforward calculation shows that E,exp(sl 2 (f)) = oo for all s > 0, whenever 
A(-B) > 0. Our next result shows that this is a quite general property of Poisson func- 
tionals. 

Corollary 3.7. Let F = Y^=i In(fn) with fn £ Ll(\ n ), n G N, and assume that 

/(®iil/nJ) CT dAl CT l < oo for all a G Il( m,...,ne) with n x , ...,n^eN 

and iaai E(^^L 1 |/ n (/n)|)^ < oo for all I G N. In addition, suppose that f n > /or a// 
n6N and t/iai taere zs an Uq > 2 wrain ||/ no ||no > 0- TTien Eexp(sF) = oo for all s > 0. 

Proof. The idea of the proof is to show that the number of partitions we sum over in the 
formulae (13.91) and (13.101) is rapidly increasing in I. For I G N with 6|£ (i.e. I a multiple 

(2) ~ 

of 6) let 11}' C n> 2 (2,...,2) be the set of partitions of [21] that can be constructed 
in the following way. First, the odd numbers in [2£] are partitioned into blocks of size 
6. Then we form £/6 — 1 blocks of size two from the even numbers of [2£] such that 
all partitions from n> 2 (2, . . . , 2) that contain the subpartition constructed so far must 
belong to Il>2(2, . . . , 2). Finally, we combine the remaining |£ + 2 even numbers of 
into blocks of size two. It follows from this construction and a short computation that 

|n ( 2 ), > I! (§1 + 2)! > l\ (1/6)! 

' 1 ' " (£/6)!(6!)l(f + 1)!2§ +1 " W 
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Note that we do not take into account here the different possibilities of forming the first 
£/6 — 1 blocks from the even numbers of [2£]. The previous inequality implies that 

|n> 2 (m,...,n,)| > |n> 2 (m,...,n,)| > |nf | > £ -j^ (3.11) 

if fix > 2, . . . , Hi > 2 and £ G N is such that 6\£. 

Each a with £ G N and Q\£ induces for every k > 2 a partition a G Il> 2 (/c, . . . , k) 

of in the following way. First, one partitions {1 + jk : j = 0, ...,£ — 1} as the odd 
numbers in a and then for each i = 2, . . . , k the sets {i + jk : j = 0, ...,£ — 1} as the even 
numbers in a. We denote the set of these partitions by IT^. 

Due to the assumptions that ||/ no || no > and that f no is non-negative and the structure 
of II^ there must be a constant c > such that 

/(OxdAH> c * (3.12) 

for all £ G N with 6|£ and a G Il^ ri0 \ It follows from Corollary 13.51 and the estimates (13.111) 
and fl3TT2j) that 

RF«>*» and -, f (F)>« C ' 
for all £ G N with 6\£. This implies that Eexp(sF) = oo for all s > 0. □ 



4 Asymptotic behaviour of moments and cumulants 

In this section we consider Poisson processes i] t with intensity measures X t : = t\, t > 0. 
We are interested in functionals of r] t that can be represented as F t — g(t) ^ f(xi, ■ ■ ■ , x m ) 
with the sum running over all m-tuples of distinct points of rj t for some integer m > 1. 
This setting is taken from [161 Section 5], where a central limit theorem for F t as t — > oo 
is derived. We generalize this to a multivariate setting and investigate the asymptotic 
behaviour of such functionals as t — > oo. More formally, fix £ > 1 and for i — 1, . . . , £ let 
mi G N, G Lj(A m< ) and & : (0, oo) ->• R such that ^ for all (or at least for all 
sufficiently large) t > 0. Now define 

F^:=gi(t) y/«(x 1 ,...,x m Jr ? ; mi) (d(a; 1 ,...,x m J), t > 0. (4.1) 

By (123]) . we have 

EF« = | /»(xi, • • .,x mi ) A^(d(x!, . . . ,x m )). 

For n = 1, . . . , mj, define 

f^(x u ...,x n ):= ( m ) f f^(x 1 ,...,x n ,y 1 ,...,y mi _ n )\ m *- n (d(y 1 ,...,y mi _ n )) (4.2) 
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and denote by I n>t the n-fold Wiener-Ito integral with respect to m. We claim that, 
P-almost surely, can be written as 

mi 

F f W = Eif } + 9i (t) J2 tmi ~ ni nAf^)- (4-3) 

n=l 

Indeed, if ffi £ L 2 (X n ) for all n < then is square-integrable and (14.31) is just a 
special case of the chaos expansion ( 12.41) . cf. Lemma 3.5 in [16]. The L 1 -version can be 
derived by approximation or by a direct calculation (just plug (14.21) into (12. 2h and observe 
that all resulting terms cancel out, except the integral representation (14.11) of F^). 

Write || • || n for the norm and (•, • ) n for the inner product in L 2 (X n ), and assume 
again that /„ £ L 2 (X n ) for n < rrii. Equations (12.31) and (14.31) imply that 

mi „ mt 

Var[F t W ] = gi {tfY,^- 2n n\ / (f®f 1 d\? = g t (t) 2 J^t 2m ^ n\ \\f^\\ 2 n (4.4) 



n=l n=l 

and that 

min{mi,mj} 

— yi\ b )yj\ L i 

n=l 



Cav[F^,F? ) ] = g i (t)g j (t) £ «! /„ W >»- (4-5) 



The variances, covariances and mixed moments and cumulants of higher order show the 
following asymptotic behaviour as t — > oo: 

Theorem 4.1. Assume that /(<8>i=i|/^ \)a dA' CT ' < oo for all a £ IT(mi, mi). Then 

»=i yn y aen =2 (i,...,i) 

and 

7 (F t (1) -EF t (1) ,...,Ff -Eif } ) 
lim ^ — ^ = 0, t > 3. 4.7 

nti 0i(*)* m4 - 1/a 



Remark 4.2. Note that the right-hand side of (14.61) vanishes for odd I. Moreover, 
7(F t (1) - EF t (1) ) = and for t = 2, the left-hand side of flU) coincides with that of 
( O]) and equals (/ x (1) , f[ 2) ) v 

Proof of Theorem \4.1\ We can assume without loss of generality that gi(i) = 1. Due 
to the special structure of F t and the kernels of its chaos expansion, the integrability 
assumptions on (8>f =1 /^ imply that J ((8>f = i I /"•< \)<r dA'' 7 ' < oo for all a £ II(ni, . . . , ni) and 
1 < rii < rrii, i — 1, . . . , I. The latter is condition (13.81) for the functions fn} in ( 14.31) and 
allows us to apply Theorem 13.11 which yields 

e n(if } - } ) = e e n t^i^fS) 

i=l l<ni<mi,...,l<ni<mi i=l 

E E [(®t=it mi - ni f$)«d\w. 

l<ni<mi,...,l<n e <m e a-eU >2 (n 1 ,...,n e ) 
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On the right-hand side, each summand has order t£( m !- n i)+M. Because of \o~\ < [(Yl n i)/2\ 
and ^2rii > £ we have XX m « ~ n i) + \°~\ <J2 m i— \(J2 n i)/2\ < — so that 

the maximal order is at most t^" 1 * - ^/ 2 !. p or even I this is obtained if and only if 
ri\ = . . . — ni — 1, and the partition a satisfies | J\ = 2 for all J G a. Exactly these sum- 
mands remain as t — > oo since they have the same order as the denominator in (14.61) ; other 
summands vanish as t — > oo. If £ is odd, the numerator has at most order t^ mi ~^ +l ^ 2 (in 
fact the order is attained) and the denominator has order t^ mi ~ l l 2 so that the expression 
vanishes in the limit. 

For the cumulant 7(F 4 {1) - EF 4 {1) , . . . , F t (e) - EF t {e) ), i > 3, we obtain by Theorem O 
the expression in the second line of (14 . 8 [) where this time the inner sum only runs over 
all partitions a G IT> 2 (ni, . . . , rii). Since IT> 2 (1, . . . , 1) D n =2 (l, . . . , 1) = for I > 3, all 
summands have a lower order than the denominator in (14.71) and vanish as t — > oo. □ 

In the next result we take F} 1 ^ = ■ ■ ■ = F^ = F t with 

F t :=g(t) J f(x 1 ,...,x m )4 m) (d(x 1 ,...,x m )), t>0, (4.9) 

as in (14. ip . where gi — ■ • • — ge — g, = ■ ■ ■ = = f and m\ — • ■ ■ — rri£ — m. Since 
n =2 (l, • • • , 1) has cardinality 

(£-1)!! := {£- l)(£-3)-... -3-1 

for even £ > 2, Theorem 14.11 implies the following result. 

Corollary 4.3. Assume that /(|/|® £ ) CT dA'°"' < oo for all a G U(m, . . . , m) and that 
|| A IK >0. Then 

E(Ft-EF t y ((£-l)l\, if £ is even, 

lim 7T- = < 

(Var[F t ])3 [0, if £ is odd, 

and 

l im ( Ft-EF t \ Jl, if £ = 2, 
^ le \y^4Ft\) [0, if £^2. 

5 Central limit theorems 

In what follows, we assume the same setting as in the previous section. More precisely, fix 
£ G N, let F t {1) , F t {l) be defined as in (jUJ) and assume for each % < £ that G L 2 (X n ) 
for n < rrii. We shall at first show how the results of the previous section lead to a 
multivariate central limit theorem via the method of moments. Let us define 

F« : = gi (t)-H-^- ll2 \F^ - EF t (i) ) (5.1) 

and note from (14.31) that 



F(W/ 2 gr^ n ,(/«). (5.2) 



n=l 
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Furthermore, by (14. 5p . we have the asymptotic covariances 

Cq:=^Cw[Fi < >,Ft l) ] = {fjf>,fP) 1 = [ fl%)fP(x)\(6x), i,j £ {1, ...,£}. 

Proposition 5.1. Let N be an £- dimensional centred Gaussian random vector with co- 
variance matrix (CV,-)ij'=i,.--/ and assume that J |(<S^ = i/^) CT | dA' CT ' < oo for all a £ 
II(mi, . . . , me). Then [F^\ . . . , F^\ converges in distribution to N . 

Proof. Observe first that ^(F t {i) ) = EF t (i) = for 1 < i < t and 

7 (F 4 (i) , F t U) ) = Cov[F 4 (i) , F t U) ] -> dj as t -> oo 

for any 1 < i, j < I. Now fix integers k > 3 and 1 < i\ < . . . < ik < £, and consider the 
joint cumulant ^(F} 11 ^, . . . , Ft). By homogeneity and (15. ip it follows that 

,fi(n) p « k) ._ 7(F™ -EF*\...,F™ -EF™) 

whence Theorem 14.11 implies that ^(F^, . . . , — > as t — > oo. The method of 
moments (or cumulants) now yields the multivariate limit theorem, cf. j2j p. 352]. In the 
univariate case the conclusion can also be directly drawn from Corollary 14.31 □ 



We now turn to a quantitative version of the multivariate central limit theorem. We 
measure the distance between two ^-dimensional random vectors X and Y by 

d 3 (X,Y) :=Bup|E0(X)-%(y)|, (5.3) 
gen 

where % is the set of all functions h £ C 3 (M^) that satisfy 

d 3 h(x) < 

dx^ dXi2 OXi^ 

Note that convergence under the (pseudo-) metric d 3 implies convergence in distribution. 
In [2], bounds are derived for the ^-distance to the multivariate normal, along with 
similar bounds using a similarly defined ^-distance. We work with the result for the da- 
distance since the covariance matrix of the Gaussian random vector is allowed to be only 
positive semz-definite (this means that some linear combinations of the components of the 
limiting random vector may be constant). A non-trivial example for such a degenerate 
situation can be found in [1]. The multivariate normal approximation of Poisson U- 
statistics in the ^-distance has been considered in [TTj . 

In contrast to the univariate results for the Wasserstein distance discussed in the 
introduction, we can derive a multivariate result only for the c^-metric since the underlying 
result in [TJ] is based on that distance. This is caused by the fact that the approaches 
used in [2] for the multivariate normal approximation, namely an interpolation technique 
and the multivariate Stein's method, require a higher degree of smoothness for the test 
functions. 

We are now ready to state the Berry-Esseen-type inequality. 



max 

l<i!<i 2 <l 



sup 



d 2 h(x) 



< i, 



max sup 
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Theorem 5.2. Let N be an £- dimensional centered Gaussian random vector with covari- 
ance matrix (Cij)ij=x g. Assume that J \fi | 3 dA < oo for every i G {1, ...,£}. T/ien 
i/iere is a constant c > suc/i £/ia£ 



ds((*l 



(i) 



Remark 5.3. For £ — 1 it is possible to replace d% in Theorem 15.21 by the classical 
Wasserstein distance dw and obtain that dw(Ft,N) < ct -1 / 2 , where N is a standard 
Gaussian random variable with variance ||/i||f and c is a constant (see [9], Theorem 7.3] 
or [THl Theorem 5.2] for a different rescaling). If ||/i||i = 0, this implies convergence in 
distribution to the constant random variable N = 0. In this situation Theorem 7.3 in [5] 
yields convergence in distribution to a higher-order Wiener-Ito integral with respect to a 
Gaussian random measure after a suitable (different) scaling. 



Remark 5.4. Theorem 15.21 also holds if F t 



(0 



are finite sums of random variables 



of type (14. ip . In fact, under some additional conditions, any Poisson functional with finite 
Wiener-Ito chaos expansion can be represented in such a way, cf. [T6] . 

We prepare the proof of Theorem 15.21 by the following lemma: 

Lemma 5.5. LetX andY be £- dimensional random vectors withKX = KY and Euclidean 
norms \ \X\\ and \ \Y\\ such that E||X|| 2 < oo, E||y|| 2 < oo. Then 

ds(X,Y) < 

{X x , 



Proof. For h G H and X 
theorem 



E||X|| 2 + E||y| 
X t ),Y ={¥,,. 



>'E\\X-Y\ 
Ye), we obtain by the mean value 



\Eh(X) - Eh(Y)\ = \E[h'(Z)(X - Y)] - E[h'(0){X -Y)]\, 

where Z = Y+U(X — Y) for some random variable U in [0, 1] and where we have used that 
the components of X — Y all have expectation zero. Applying the mean value theorem 
again as well as the Cauchy-Schwarz inequality yields 



\Eh{X) -Eh(Y)\ 



i=i ^ ' 



d 2 h(Z 



E ' ' 

i=i j=i j l 



■Zj(Xi - Yi) 



< 



ditjdui 
j=i ■> 



E X -Y 



with random vectors = UiZ and random variables Ui E [0, 1], % — 1, 
fact that h G % and the Cauchy-Schwarz inequality, it follows that 

2 



I. By the 



i=\ \j=i 



d 2 h 

du^du. 



{Z {i) )Zj < £ 2 E\\Z\\ 2 < i 2 (E\\X\\ 2 + E||F|| 2 ) 



which completes the argument. 



□ 
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Proof of Theorem \5.2[ For i £ {!,...,£} define 



and note that Cov[F^,F^] = Cij. Therefore we obtain from [TH Corollary 4.3] that 

d,((F«, ...,FW ),N) <ji2 r ' 3/2 I \fi l \x)\ 3 Xt(dx)=c 2 r^ (5.4) 

for some C2 > 0. Lemma [5.51 implies 

d,({FP,...^UFP } ...^))<A /2 Bl'\ (5.5) 

where 

A t : = ^E(F«) 2 + ^E(F«) 2 , B t := £l(F« - F«) 2 . 

i=l i=l j=l 

The first factor t4 4 is bounded in t. For the second factor we use (15. 2p to obtain that 

i=l n=2 j=l n=2 

so that A^J 2 B^ 2 < c^t~ x l 2 ^ t > 1, for some C3 > 0. Using this estimate in (15. 5p and 
combining with f ]5 .41) and the triangle inequality for 0I3, we obtain the result. □ 

Remark 5.6. The proofs of Theorem 15.21 and the univariate bound discussed in Remark 
I5.3l depend on general Berry-Esseen type inequalities for Poisson functionals from [121 [14] . 
that are proven in a slightly more restrictive setting, namely that (X, X) is a Borel space 
and /1 is non-atomic. But they are still valid without these assumptions since the proofs 
only make use of properties of the Malliavin operators that also hold in our more general 
setting as shown in [TU] . 



6 Poisson flat processes 

In this section we assume that r] t is a stationary Poisson process of k-flats (k- dimensional 
affine subspaces) in M. d , where d > 1 and k £ {0, . . . , d — 1}. This is a Poisson process 
on the space A(d, k) of all /c-flats, whose distribution is invariant under translation of the 
flats. Its distribution is determined by the intensity t > and the directional distribution 
Q, a probability measure on the space G(d, k) of all A;- dimensional linear subspaces of M. d . 
In fact, the intensity measure A t of r\ t equals 

\{-) = t\ [ l{E + xe -}U d ~ k (dx)Q(dE), (6.1) 

JG(d,k) JEJ- 

where 'H d ~ k denotes (d— fc)-dimensional Hausdorff measure, and Q is a probability measure 
on G(d,k). We let A := Ai. If Q is the uniform distribution (Haar measure), then r] t is 
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isotropic, that is, distributionally invariant under rotations. For further details on Poisson 
flat processes we refer to [T7] . 

The intersection process of order m G N is given as the set of all intersections E\ D 
• • • H E m of m pairwise different flats in r] t . To introduce our geometric functionals of the 
latter process we let C d denote the system of all compact subsets of R d , equipped with 
the Borel cr-field induced by the Fell topology, see e.g. [TF1 Chapter 12]. We consider a 
measurable family Cq C C d of sets containing the empty set and with the property that 
rBHE eC$ for all B G C$, all r > 0, and all affine subspaces E C R d . We assume that 
ip : C d — > R is a measurable function with t/>(0) =0 satisfying 

J mB n E 1 n ■ ■ ■ n E m )\ 3 \ m (d(E u . . . ,E m )) <c B (6.2) 

for all B G C$, where C B > is a constant only depending on B. (By [TT] Theorem 12.2.6] 
the mapping (Ei, . . . , E m ) ^ B C\ E\C\ ■ ■ ■ C\ E m is measurable.) We note here that ( 16. II) 
implies that A is locally finite in the sense that A({i? G k) : B C\ E ^ 0}) < oo for all 
B G C d . Since ^(0) = 0, assumption (EJ) implies the integrability of ^(BnEiH- ■ -n£ m )| p 
w.r.t. A m for any p G (0, 3]. This is enough to settle all integrability issues in this section. 
Clearly (16. 2 j) is implied by the stronger condition 

WBnEin-nE m )| < c B , A m -a.e. (£ 1; . . . ,£ m ), 5 G C d (6.3) 

for some cb > depending on B. In particular, (I6.3P is satisfied in our examples below. 
Define a random field ( t := {( t (B) : B G Cq} by 

Ct(£):=^j J ip(BnE 1 n---nE m )ril m) (d(E 1 ,...,E m )), BeC d . 

Since 7^ has only atoms of size one (by (16. ip ) we can identify rj t with its support, and inte- 
gration with respect to % corresponds to summation over all m-tuples (E±, . . . , E m ) G 
rj™ with pairwise different entries. For A,B G Cq define 

C(A,B) : = * y ( y ^in^n^n-.nEj A™- 1 (d(3>,..., 

x ( | ^(Bni? 1 n^n--n£jA m - 1 (d(^... 1 £:)))A(di? 1 ). 

If m = 1, this has to be read as 

C(A 5) = J iP(A n n E{) \(dEi). 

It can be checked directly that C(-, •) is positive semidefinite. Therefore we can consider 
a centred Gaussian field £ := {£(-8) : B G Cq} with this covariance function. 
Define 

Q(B) := r^- l ' 2 \Ct{B) - ECt(B)), t > 0, B £ C d . 
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Theorem 6.1. Let £ > 1 and B 1 ,...,B e e Cjf. Then 

d 3 ((Ct(Si), • • • , Ct(Be)), (t(Bi), Z(B t ))) < 0(5!, . . . , B t )rV 2 , t>l, 
for some constant c(B\, . . . , Bp). In particular 

{( t (B) : B £ C d } -A {£(5) : B £ C d Q } as t oo 
in t/ie sense of convergence of finite- dimensional distributions. 

Proof. This is a direct consequence of Theorem 15.21 □ 

Alternatively one can approach the central limit problem in another but closely related 
setting. Instead of increasing the intensity parameter t, we can also fix t (for simplicity we 
take t — 1) and increase the size r of the observation window. If we assume additionally 
that the considered function ip is homogeneous of degree aeR, that is 

i/j(rB) = r a i)(B), B £ C d , r > 0, (6.5) 

both approaches are equivalent. Define a random field ( r := {( r (B) : B £ Cq} with 
( r (B) = r -(m-i/2)(d-*)-«(^( r s) _ ECi(rS)). 

Corollary 6.2. ^sswme (KTBjL let £ > 1 and B 1 ,...,B e G Cjf. Then 

d 3 ((( r (Bi), Sr(B t )), (£(£i), . . . , e(B«)) < 0(5!, . . . , 5,)r-( rf - fc )/ 2 , r > 1, 
/or some constant c(B\, . . . , S^). In particular 

{( r (B) : 5 £ C rf } -A {£(£) : B £ C d } as r ^ oo 
m t/ie sense of finite- dimensional distributions. 

Proof. The special structure ( 16. ip of the intensity measure A implies the well-known scal- 
ing property 

P(7ft £ •) = P(t~ 1/(d ~% £ •), t> 0, 

where ar/i := {aE : 2£ £ rji} for a > 0. Since ^ is homogeneous we obtain for all B £ Cq 
and r > that 

&(rfl) = ^r Q / ^nr^n-nrXj^td^i,..,^)) 

ml J 

= [ 4>(B n 5i n ■ • ■ n £ m ) (d(£ l5 . . . , £ m )) , 

771. J 

where the second identity holds in distribution jointly in B. Hence, we can apply Theorem 
EZSwith gi (t) = ■■■ = ge (t) ■= {m\)- 1 t a{d - k) and then replace t by r 1/(d ~ fc) . □ 
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Remark 6.3. It follows from (14. 5p (with g^t) = gj(t) = 1/m!) that 



Cov[C t (A)Xt(B)] = Yl 



J n\((m — n)\) 



; V t (A,B,n), A,BeC d , 



where 



V t (A, B, n) := t 2m ' n J J ip(A n £1 n ■ ■ ■ n n E n+1 n-n£ m ) \ m ~ n (d(E n+1 , . . . , E m )) 
J ^{B n E 1 n • • • n E n n E n+1 n • • • n E m ) \ m - n (d(E n+1 , E m j) 



\ n {d{E u ...,E n )). 



In accordance with Theorem 16.11 we therefore obtain 

hm t^ 2m ^ Cov[Ct(A), Ct(B)} = C(A, B). 



t— >oo 



We now present a couple of examples to which Theorem 16.11 as well as Corollary 16.21 
can be applied. 

Example 6.4. Assume that m(d — k) < d. Assume further that Cq = C d and that ip is the 
(d— m(d— A;))-dimensional Hausdorff measure on ~R d restricted to C d . Then ( 16. 5 p holds with 
a = d — m(d — k). Assumption (16. 3p holds because for A m -a.e. (E±, . . . , E m ) e A(d, k) m 
the intersection E\ n • • • fl E m is either empty or has dimension d — m(d — k). This follows 
(recursively) from the argument given in [T7J p. 130]. 



Example 6.5. Assume that C d 
a = while (16. 3p holds with eg 



C d and that if>(B) = 1{B ^ 0}. Then fl§3]) holds with 
1. 



Examples 16.41 and 16.51 have been studied in [U [S] in the case k = d — 1. Our results 
add to Theorem 3.1 and Theorem 4.1 in [1] in several ways. While the latter results are 
multivariate central limit theorems for the d possible values of the number m of intersec- 
tions but a fixed (convex) test set B, we fix m but study (t(B) (respectively (i(rB)) as a 
function of B. Further we consider processes of flats and not only hyperplanes. Moreover 
we obtain Berry-Esseen-type bounds on the distance d% and can allow for a considerably 
larger class of functionals ip. It is also possible to apply Theorem 15.21 to the vector- valued 
processes arising by varying m. This would constitute a complete generalization of [I]. 
In order to avoid heavy notation we have refrained from doing so. 

We continue with further examples of functionals i/j satisfying (16. 3p and (16.51) . The 
convex ring lZ d is the system of all (possible empty) unions of convex and compact subsets 
of R d . 



Example 6.6. Assume that Cq = TZ d and that ip is the intrinsic volume V a , where 
a G {0, . . . , d}, see e.g. [T7|. Then (16.51) holds. Assumption (16. 3p follows from the fact 
that V a (B n E) < V a (B) for any convex and compact B C M. d and any affine subspace 
E C M. d . By additivity of V a (see e.g. pH Section 14.2]) the inequality (Q can be 
established for the whole convex ring. 



In contrast to the previous examples, the next functionals are not additive. 
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Example 6.7. Assume that Cq = TZ d and a G {0, . . . , d — 1}. Let Q (A, ■) be the support 
measure of A G IZ d , see [TTl Section 14.2]. This is a signed measure on the product of 
W 1 and the unit sphere S d_1 such that Q a (A, W 1 x § d_1 ) = V Q (A). Fix a measurable set 
U C and assume that 



where N(A) is the unit normal bundle of A. This consists of all pairs (x, u) G M. d x S^ -1 
that occur as unique nearest point and associated direction of a point in the complement 
of A, see [6]. The homogeneity (16. 5p follows from the homogeneity of the (non-negative) 
measure l{(x, u) G N(A)}Q a (A, d(x, u)), see [6j Proposition 4.9]. Assumption (16. 3p fol- 
lows similarly as in Example 16.61 from the additivity of Q a (A, •) in A G 7?. d . 

Example 6.8. Consider the case where Cq is the space of compact convex subsets of M. d , 
fix a G {0, . . . , d} and /3 > and let -0 be V£, the power /3 of the intrinsic volume of order 
a. In the case a = 1 and /3 = n G N, ^ corresponds to the n-th chord-power integral, 
which is frequently studied in integral geometry, cf. [T7J Chapter 8.6]. Clearly, (16. 5 p is 
satisfied with aj3 there and assumption (16.31) follows as in Example 16.61 from the fact that 
V^(BdE) < V£ (B) for any convex and compact BcM d and any affine subspace E C M. d . 

Remark 6.9. If ij) > 0, then Cov[Ct(A), ( t (B)] > and C(A,B) > for all A, B G C^. 
This is the case in Examples I6.4[ I6.5[ and 16.81 Taking as Cq the system of convex sets, 
this is also the case in Examples 16.61 and 16.71 If additionally m > 2, Corollary 13.71 shows 
that the moment generating functions of the functionals under consideration do not exist. 
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